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I. INTRODUCTION 

Quantum correlations between particles constitute an 
important topic in quantum information research. Uni- 
versal schemes for quantum computing allow the forma- 
tion of arbitrarily correlated states from simple product 
states of two-level quantum systems, and via sequential 
application of one- and two-particle interactions, a wide 
variety of many-body interactions may be effectively sim- 
ulated [ll-0]- Particular interest has been devoted to 
topological matter which shows correlations but no lo- 
cal order parameter [1-0|- The Laughlin wave-function 
describing the fractional quantum Hall effect is an ex- 
ample of a topological state of correlated electrons, while 
high-Tc superconductivity has been proposed to originate 
from the resonating valence bond (RVB) state, which is 
a superposition of different spin singlet configurations. 

While demonstration of topological matter in natural 
physical materials is so far restricted to the quantum Hall 
effect, a number of proposals exist to synthesize states of 
matter, which display the desired correlations. One may 
distinguish two kinds of implementation: i) synthesis of 
states of many particles which are formally equivalent to 
topological matter candidates, ii) synthesis of effective 
Hamiltonian interactions, which have the desired topo- 
logical states as ground states or eigenstates. The first 
kind of implementation may be used to verify by mea- 
surement the value of certain correlations and effects of 
certain operations on the system, while the second kind 
may be used to illustrate phase transition dynamics, dy- 
namical response, and robustness to external perturba- 
tions. 

Topological states on a two dimensional array of spin 
1/2 particles or qubits, are interesting candidates for 
quantum memories stored in collective states of many 
particles. Quantum information is stored in states with 
large topological degeneracy, but with an excitation en- 
ergy gap which protects against local disturbances d-Q . 
Even without interactions, specifically ordered states of- 
fer error protection, if one is able to probe the multi- 
particle constraints and thus identify single particle er- 
rors when they occur in quantum computing [9| and com- 
munication [10| and suitably correct them. 

Periodic arrangements of atoms in optical lattices 



physically display the superfluid-insulator dynamics of 
the HublDard Hamiltonian [ill, 113) ^"^^ lattices with a 
single trapped atom per lattice site have been proposed 
as simulators of Ising and Heisenberg interacting spin 
models [l^, and more elaborate lattice gauge theories 
[3 [l^, where the geometry of the potential wells in 
an optical lattice are here in immediate correspondence 
with the mathematical models of particles or spins on a 
square, triangular or honeycomb lattice. 

In fTG*!, Paredes and Bloch offer an excellent review of 
the properties of a variety of topological states together 
with proposals for implementation of these states with 
neutral atoms trapped in optical lattice potentials. In 
particular, they focus on a single plaquette, i.e., an ar- 
rangement of four spin 1/2 particles at the corners of a 
square as illustrated in Fig. [Tfa). If one represents the 
spin up (down) states by filled (empty) sites of hard core 
bosons, one obtains the state illustrated in Fig. (rfb). 

Tunnelling among wells combined with controlled 
ground state collisions [ll] and long range dipole interac- 
tion between atoms excited to high lying Rydberg states 
[l3 [ have been proposed to engineer both unitary and 
dissipative evolution leading to formation of topological 
states in the lattice system. 

In this work, we will show how Rydberg blockade in 
small atomic ensembles in conjunction with qubit encod- 
ing in the collective atomic population in different in- 
ternal atomic states ^17] can be used to simulate simple 
instances of topological matter. In particular, we will 
show that a few hundred atoms with five stable internal 
states can be used to simulate Laughlin-like states, RVB- 
states, and flux and charge quasi-particle excited states 
of a simple plaquette system. Such states provide an in- 
teresting stepping-stone for experimental demonstrations 
and investigations of topological effects. The collective 
qubit encoding in atomic ensembles is formally equivalent 
to photonic quantum computing schemes. Multi-photon 
states have been produced with correlations mimicking 
the spi n correlations of simple models of topological mat- 
ter [isl Il9| . The photonic implementations use squeezed 
light sources, and detector coincidence signals verify the 
values of quantum correlations in a heralded (and de- 
stroyed) component of the light. In contrast our atomic 
scheme generates the states on demand and permits se- 
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FIG. 1. (Color online) (a) A plaquette of four spins in an 
optical lattice, (b) Equivalent plaquette of hard core bosons 
with a filled (empty) site representing the spin up (down) 
states in part (a) of the figure. 

quences of subsequent measurement and interaction op- 
erations to be carried out on the system. 

In Secini we discuss the main features of the plaquette 
model that we will implement and study in this paper. In 
Sec. llin we review the Rydberg blockade mechanism and 
we discuss and analyze the schemes for collective qubit 
encoding and manipulation needed to simulate topolog- 
ical states in small atomic ensembles. In Sec. IIVI we 
present effective protocols to engineer particular topolog- 
ical quantum states in collective atomic degrees of free- 
dom, and we perform a variety of operations on these 
states. Section IVl concludes the paper. 

II. THE PLAQUETTE 

A simple instance of topological matter can be simu- 
lated by four spin 1/2 particles arranged in the corners of 
a square and enumerated 1-4. Figure [IJa) shows a clas- 
sical spin configuration with spins 1 and 3 in spin down 
states and spins 2 and 4 in spin up states. Spin singlet 
states of any short range pair of spins, i.e., spins on the 
same side of the plaquette such as spins 1 and 2, can be 
defined, and the resonating valence bond (RVB) states 
are superpositions of states where the spin singlets occur 
for the pairs (1,2) and (3,4) and for the pairs (1,4) and 
(2,3). These RVB states are examples of topological spin 
liquids. 

In Fig. mb), the spin down particle is replaced by an 
empty site, and the spin up particle is replaced by a site 
occupied by exactly one particle. The Hilbert space of 
four spins with a vanishing total z-component is thus 
equivalent to the space of two hard core bosons dis- 
tributed among four spatial locations with at most one 
boson per site. The state shown is thus in one-to-one 
correspondence with the state shown in part (a) of the 
figure, but the quantum degree of freedom is here associ- 
ated with the occupation of the sites, i.e., by the spatial 
location of physical particles, rather than by the internal 



spin degree of freedom of space fixed particles. Note that 
it is a natural property of this representation that a sin- 
gle particle occupying an odd superposition state of two 
lattice locations is equivalent to an entangled spin singlet 
state of the corresponding two spin 1/2 particles. 

In [l^, the analogy of the hard core boson occupa- 
tion degree of freedom and the spin 1/2 internal states is 
discussed in detail, and it is proposed how laser double 
well potentials along two spatial directions can be used 
to define the plaquette and how tuning of the on-site in- 
teraction and the tunneling of particles may be used to 
engineer at will a variety of different states in both the 
spin and in the hard core boson representation. It is for 
example possible to allow tunneling with amplitude tx 
along the 1-2 and 3-4 bonds, causing a super exchange 
interaction Hs = Jx{Pi2 + ^34), where Pij is the pro- 
jection operator on the spin singlet state of spins i and 
j, Jx = 4:t'^/U with U the on-site interaction energy. 
The product state of (1, 2) and (3, 4) spin singlets is then 
gapped from all other states of the four spins, and by 
lowering the potential barrier in the orthogonal direction 
and thus increasing the tunneling amplitude, ty, along 
the 1-4 and 2-3 bonds, one provides a Hamiltonian 

Hs = Jx{Pl2+p3i) + Jy{Pli + P23), (1) 

which adiabatically connects the product state of two 
simple spin singlets with one of the RVB states. In the 
hard core boson representation (removing the spin down 
particles) the RVB state of the remaining two spin up 
particles is also equivalent to the Laughlin state, and 
in 16] it is proposed how further use of external inho- 
mogeneous perturbations can engineer quasi-particle ex- 
cited states and simulate their anyonic behavior associ- 
ated with the excitations being propagated around all 
corners of the plaquette. 

III. RYDBERG BLOCKADE AND COLLECTIVE 
ENCODING OF QUBITS 

Instead of coding memory bits into separate particles, 
it has been proposed to approach scalability of quantum 
information processing by making use of the rich inter- 
nal level structure found in many atomic and molecu- 
lar quantum systems [13, HO]- In these proposals it is 
suggested to use ensembles of identical particles and to 
encode quantum information in symmetric states of the 
ensemble with different values of the collective popula- 
tion of the different internal atomic or molecular states. 
To be more specific, an arbitrary A^-bit register state 
\b1b2 ■ ■ - bpf), with hi = 0,1, is encoded in the collec- 
tive symmetric state of K identical quantum systems 
with (N + 1) internal levels {\i) ,i — 0, . . . , N} if one lets 
precisely hi particles populate the internal state |i), see 

Fig. [21 while K—J2i=i bi particles remain in the so-called 
reservoir state |0). Values of hi that are larger than unity 
are not compatible with the binary logic of qubits, but 
access to an extra quantum level associated with each 
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FIG. 2. (Color online) Collective encoding of a four qubit reg- 
ister with zero or unity collective occupation in four different 
internal atomic states. The state shown is equivalent to the 
hard core boson plaquette illustrated in Fig. [Hb). 



qubit can lead to very efficient operations, as we shall 
see in the following. 

The system may be initialized with all atoms prepared 
in the state |0), e.g., by optical pumping, and provided 
the lasers and other means used to excite and control the 
ensemble act identically on all atoms, the state will re- 
tain its permutation symmetry for all times. Note that 
for atoms at rest, e.g., trapped in an optical lattice po- 
tential, phase factors e*'^ '"^ arising from travelling wave 
laser fields which excite the atoms at their different spa- 
tial locations can be absorbed in a redefinition of the 
atomic internal state coherences. To force the system to 
remain in the subspace of meaningful quantum register 
states we propose to apply the Rydberg blockade mech- 
anism ^li, i22i] , which uses the strong long-range interac- 
tion between highly excited atoms and the associated en- 
ergy shifts to prevent atoms in the vicinity of an already 
excited neighbor from being transferred into a Rydberg 
state. For a review of quantum information processing 
ideas with Rydberg excited atoms and relevant physical 
parameters, see [23[. 

In an ensemble of atoms, which are all within reach of 
these strong interactions, extending out to 10 /im separa- 
tion, a resonant laser process towards a definite Rydberg 
state may produce an excited superposition state 

1 ^ 

1^) = 7f E . . . 0,-ir,0,+i . . . Oa') (2) 

in which all atoms are excited with the same probabil- 
ity amplitude, but the laser field resonant with a single 
Rydberg excitation is detuned for states with two or more 
Rydberg excited atoms because of the interaction energy 
between excited atoms. Taking advantage of this mecha- 
nism one can drive precisely one atom from the initially 
macroscopically populated state |0), via a Rydberg state 
|r), into an initially empty internal state, and thus pre- 
pare a collective state with precisely one atom or no atom 
in any of the register states 

The coherent excitation towards the Rydberg state 
may be terminated so that a superposition of a state 
with zero and one Rydberg atom is produced, and in 
this way, we can prepare quantum superposition states 
of the qubit register. The excitation to a Rydberg level 
from one internal state can also prevent excitation from 



another internal state, and in this way two-qubit entan- 
gling operations can be implemented through the same 
mechanism In fact, the universal power of the collec- 
tive encoding quantum computing allows us to synthesize 
any state and to implement any unitary transformation 
within the atomic ensemble. 

It is important to point out, that the operations men- 
tioned are carried out by collective addressing, and the 
exciting laser pulses only need to be adjusted to the tran- 
sition frequencies, the coupling strengths and the polar- 
ization selection rules of the atomic transitions. A practi- 
cal implementation in alkali atoms may thus apply a ho- 
mogeneous magnetic field, which splits the different Zee- 
man sublevels of the atomic hyperfine ground state struc- 
ture, and thus the excitation towards the Rydberg states 
and lower lying optically excited states can be made state 
specific via the transition frequencies (T7| . 



A. Entanglement generation with and without 
interactions 

We noted already that the universal quantum com- 
puting capacity of the collective coding scheme and the 
Rydberg blockade mechanism permit the generation of 
any state of the system. In this subsection we wish to 
point out that some states are actually more easily syn- 
thesized in an alternative manner than by processing a 
sequence of binary logic quantum gates. As mentioned 
in Sec. |lTl a single hard core boson in an odd super- 
position state of two spatial locations is equivalent to a 
spin singlet state, i.e., an entangled state of two local- 
ized spin 1/2 particles. With our collective encoding in 
Fig. [21 we similarly observe that a logical state like |1000) 
with one atom in state 1, and no atoms in the other reg- 
ister states, can be transformed into any superposition 
cillOOO) + C2IOIOO) + C3IOOIO) -I- C4IOOOI) by simply ap- 
plying (Raman) laser fields onto the entire ensemble that 
transfer non-interacting atoms from the internal state |1) 
into the superposition Ci\i). 

It may seem puzzling that without use of any interac- 
tion between the atoms, it is possible to turn a classi- 
cal register state into a quantum state corresponding to 
entangled qubits. In registers with individual quantum 
systems encoding of the different qubits such operations 
require direct or indirect interactions between the par- 
ticles. But we recall that the "classical" register state 
1 1000) is indeed a multi-particle entangled W-state of the 
atomic ensemble, where every atom has the same ampli- 
tude to populate the internal state |1). We observe that 
this atom-atom entanglement can be freely transformed 
into qubit-qubit entanglement within the register sub- 
space of only a single bit value of unity. While we need 
the Rydberg blockade to manipulate the collective pop- 
ulation between zero and unity in the register states, in 
the conventional encoding, one-bit unitary gates do not 
need any interaction. 

For our purpose, it follows from the above that pairs 
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of spin singlets of spins (1,2) and (3,4) are readily ob- 
tained from the state shown in Fig. [5] by simple Ra- 
man processes, while the superposition of (1,2), (3, 4) and 
(1, 4), (2, 3) singlet pairs in the RVB states could be non- 
trivial two-particle states and hence deserve special at- 
tention. We will now present a general criterion for which 
two-atom states can be transformed into each other by 
simple Raman processes, i.e., by simple linear transfor- 
mations on the single atom level structure. 

Symmetric collective states of if ^ 1 atoms with two 
atoms in the register states may be conveniently repre- 
sented as 



N N 

EE 

2 = 1 j=l 



ajc, 



.,at|00. 



..0), 



(3) 



where the a| operators are bosonic creation operators of 
an atom in level i. The bosonic property does not rely 
on the use of bosonic atoms in the ensemble, but merely 
reflects the permutation symmetry among the atoms, c 
is a matrix of coefficients determining the state, and for 
uniqueness, we require Cy = Cji. A Raman transition be- 
tween levels i and j (where i and j may refer to the same 
level) is described by a single atom Hamiltonian Hi = 
^(5b)(*l which gives rise to the multi-atom, 

second quantized Hamiltonian H — h{ga}jai + g*a\aj). 
In the Heisenberg picture, this Hamiltonian causes a lin- 
ear transformation among the bosonic creation operators, 
and by choosing suitable combinations of Raman transi- 
tions, it is possible to implement any linear transforma- 
tion of the form d\ — UijOj, where [/ is a unitary 
matrix (the unitarity preserves the bosonic commutator 
relations [ai,dj] = 6ij). This transformation transforms 
the ensemble state li/ji) into 



N N N N 



1^/) = E E E E 4Uuc,jU,kai\00 . . . 0). 



(4) 



i=l j=l k=l 1=1 



We can thus transform a symmetric ensemble state {ipi) 
with expansion coefficients c into a state with coefficients 
c without interatomic interactions if and only if we can 
find a unitary matrix U so that 



c = cU, 



(5) 



where T denotes the transpose. 

According to the Takagi decomposition [2^, it is al- 
ways possible to write an n x n complex symmetric matrix 
A a.s A — VYiV'^ , where ^ is a unitary matrix and E is a 
real nonnegative diagonal matrix. The diagonal elements 
of E are the positive square roots of the eigenvalues of 
AA* and the column vectors of V are a set of orthonor- 
mal eigenvectors of AA*. It thus follows that two states 
with coefficients c and c can be transformed into each 
other by use of independent atom (Raman) transitions 
if and only if the Hermitian matrices cc* and cc* have 
identical eigenvalues. (For c = V'SV'^ , c — V'SV'^ , and 
S = E, the solution of © takes the form U = V*V'^.) 



This simple analysis is useful to assess which gate im- 
plementation of quantum algorithms can be most easily 
carried out in the collective coding scheme, e.g., on re- 
mote, non-interacting particles. In the next section we 
will apply the criterion ([5]) and find simple means to pro- 
duce different plaquette states. 



IV. PREPARATION AND MANIPULATION OF 
TOPOLOGICAL STATES 



Hard core boson equivalence of resonance 
valence bond states 



As specific examples, let us assume N = A and consider 
the hard core boson states equivalent with spin singlet su- 
perposition states by the identification of the spin singlet 



creation operator s. 



|a>_) 



('Sl2'S43 ■ 



''14*23 



)|0000) = (|1001) + 10110) 

- 10011) - |1100))/2 (6) 



and 



{4,sl, + sl43)\0m/V3 - (11001) + 10110) 
|1100) + 10011) - 2|0101) - 2|1010))/VT2. (7) 



These are the resonating valence bond states of the sim- 
ple plaquette described in ^] . The coefficient matrix for 
|$_) is determined as 



c = — 
4 



0-101 
-10 10 
10-1 
10-10 



(8) 



and the eigenvalues of cc* read 0, 0, 1/4 and 1/4. |-0i) — 
1 1100) is a simple state to prepare by transfer of two 
atoms via the Rydberg state to states |1) and |2), and 
with the coefficient matrix 



10 

10 







(9) 



we find that cc* has the same eigenvalues, 0, 0, 1/4 and 
1/4. This implies that the state |<I>_) can be prepared 
from = 1 1100) by a simple transformation, and we 
compute the corresponding unitary single atom evolution 
operator. 



U 



1 

71 



-10 10 
10-1 
10 10 
10 1 



(10) 



This matrix represents a 7r/2 pulse between levels 1 and 
3 and a 7r/2 pulse between levels 2 and 4 with suitably 
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chosen phases, which appears because l^*-) is, indeed, a 
direct product of a singlet between spins 1 and 3 and a 
singlet between spins 2 and 4 [l6l |. 

We next consider the state |$+), which has the coeffi- 
cient matrix 



4\/3 



f -2 f 

1 f -2 
-2 f f 

f -2 f 



(11) 



and the eigenvalues of cc* read 0, 1/12, 1/12, and 1/3. 
For "classical" register states with only one term, the 
eigenvalues are always either 0, 0, 1/4, and 1/4 or 0, 0, 
0, and 1, and it is thus impossible to prepare |<i>+) from 
such simple states without using the Rydberg blockade 
(or another interaction between the atoms). We can still, 
however, use our mapping criterion, and, e.g., verify that 
the state a/2cii|2000) + (c23 + C32)|0110) (where |2000) is 
the state with two atoms in level |1)) leads to the eigen- 
values 0, |ciip, |c23p, and |c32p. Choosing cn = l/\/3 
and C23 = C32 = l/(2\/3), we thus find the same eigenval- 
ues as for \^+). This implies, that we can find the single 
particle unitary transformation by which these two states 
are converted into each other, and it reads 



U ■ 





I 


—I 




—I 


1 


—i 


1 


i 


-1 


2 


i 


1 


—i 


-1 




1 


1 


1 


1 



(12) 



The state a/273|2000) + ^17310110) can be prepared 
as illustrated in Fig. [S] In part (c) of the figure, we use a 
composite pulse sequence consisting of three pulses with 
constant amplitude and phase, which was proposed in 
[23 |. The ability to control the transfer of the atom in 
the Rydberg level into |1) dependent on whether |1) is 
initially populated by one or zero atoms ((c) to (d) in 
the figure) stems from the fact that the transition be- 
tween a ground state level and a Rydberg level, due to 
the blockade effect, acts as a two-state transition with a 
Rabi frequency proportional to the square root of the to- 
tal number of atoms occupying the two levels. The Rabi 
frequency on the |1) to |r) transition is thus a factor of \/2 
larger in the former case than in the latter case. We also 
investigated whether the transformation achieved with 
the composite pulse sequence could be achieved with only 
two pulses, but restricting the parameters so that the 
transfer is perfect for the state with no atoms in |1) ini- 
tially, we obtained a numerical maximum of the trans- 
fer fidelity for the state with one atom in |1) initially of 
0.7337. The three pulse scheme is hence necessary. 

We note here the convenience of using an intermedi- 
ate state, which has double occupancy of a single atomic 
level and hence does not belong to the class of hard core 
boson states. We imagine, that further exploration of 
the ladders of occupation number eigenstates may lead 
to a variety of schemes for efficient state synthesis and 
manipulation. 



B. Braiding 

Excitations of topological states correspond to anyonic 
quasiparticles, and we next describe a method to demon- 
strate the phase factor that arises when two such anyons 
are created and one is moved around the other one. Note 
that our qubits are not localized in space, but we may 
nevertheless regard the process as braiding due to the 
analogy between our system and a plaquette in an opti- 
cal lattice. As in nm^ we consider the state 



|n) = i=(|ioio) + 10101)), 



(13) 



which is a minimum instance of a string-net condensate. 
If two Rydberg levels \r) and \r') are available and the 
presence of a single atom in \r) shifts the energy of |r') 
significantly, this state can be prepared from |1010) as 
follows. We first apply a Raman pulse between |1) and |2) 
to obtain the state (|1010) + |0110))/V2. We then apply 
(i) a TT pulse between |1) and |r), (ii) a tt pulse between 
|3) and |r'), (iii) a tt pulse between |4) and |r'), and (iv) 
a TT pulse between |1) and |r). This sequence moves the 
atom in |3) to the level |4) if |1) is unoccupied and hence 
produces the desired state. With only one Rydberg level 
available, we can use the procedure in Fig. [3] to prepare 
the state (|2000) -I- |0101))/V^ and then apply a tt/V2 
pulse between |1) and |r) followed by a tt pulse between 
|3) and \r) to transfer a single atom from |1) to |3) in 
the first term. The reasoning behind this procedure is 
that we need to include a pulse sequence, which utilizes 
the dependence of the Rabi frequency for a transition 
between a register state and the Rydberg level on the 
total number of atoms occupying the two levels in order 
to obtain a qubit controlled nonlinear effect with only one 
Rydberg level. On the other hand, we can not transform, 
for instance, |2000) into |1010) or vice versa only by use 
of Raman transitions, and we thus use the Rydberg level 
a second time to achieve the desired state. 

If we denote the Pauli spin operators of the ith qubit 
by (erf , erf , erf ), we obtain a pair of so-called chargelike 
quasiparticle excitations if we apply erf, erf, erf, or cr| to 
!□). A pair of fiuxlike quasiparticle excitations, on the 
other hand, is obtained if we apply erf , erf , er| , or er| to 
!□) [l6|. These excitations have anyonic character, i.e., 
they accumulate non-trivial phase factors when moved 
around each other. In [l6| , it is thus proposed to move a 
chargelike quasiparticle excitation around the plaquette 
by subsequent application of er^ operators and study the 
dependence of the outcome on the prior excitation of a 
pair of fiuxlike quasiparticles, of which one is located at 
the center of the plaquette. 

Specifically, application of the sequence of operators 
erferf erf erf creates a pair of chargelike quasiparticle exci- 
tations, moves one of them around the plaquette, and an- 
nihilates the pair again, which altogether transforms [□) 
into !□). On the other hand, if we first apply erf to cre- 
ate a pair of fluxlike quasiparticle excitations erf !□) and 
then apply the cr|cr|er|erf sequence, we obtain — erfjD) 
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FIG. 3. Preparation of the state -\/2/3|2000) + \/l/3|0110) by use of a single Rydberg level. Full arrows are vr pulses, dashed 
arrows are Raman pulses, doublearrows are composite pulses, and |e) is an excited level, (a) Starting from the register state 
jOOOO), we apply two tt pulses in the indicated order to move precisely one atom from |0) to |1) as explained in Sec. IIIII (b) We 
apply a Raman pulse to obtain the state 1^2/3] 1000) + \/l/3j0100). (c) We first transfer one atom from |0) to |r) and then 
apply a composite pulse between levels |1) and \r), which transforms a state with one atom in |1} and one atom in \r) into a 
state with two atoms in |1) and no atoms in |r) and leaves a state with no atoms in |1) and one atom in |r) unchanged (see 
text), (d) We apply a n pulse to interchange the populations in \r) and |3), which leads to the desired state (e). 
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the second is the state of the four register qubits. (ii) 
Apply erf if and only if the control qubit is in the state 
|0) to obtain 

|0) (g) (|1010) - |0101))/2+ 11) (g) (11010) + 10101))/2. (15) 

(iii) Apply the operator afa^ (T2(yf , which changes the 
sign of the first term but not of the second, i.e., 

-|0)(8)(|1010)-|0101))/2+|l)®(|1010) + |0101))/2. (16) 

(iv) Apply af if and only if the control qubit is in the 
state |0) to obtain 



(-|0) + |l))(8)(|1010) + |0101))/2. 



(17) 



FIG. 4. Braiding of a chargelike quasiparticle around a flux- 
like quasiparticle. (a) The plaquette seen from above, (b) Ap- 
plication of af creates a pair of fluxlike quasiparticles (black 
circles), (c) af creates a pair of chargelike quasiparticles 
(white circles), (d-e) and erf move one of the chargelike 
quasiparticles. (f ) af annihilates the two chargelike quasipar- 
ticles. Finally, the two fluxlike quasiparticles are annihilated 
by applying af, and we return to (a) except for a phase change 
of n. 

with the desired change in sign. This is illustrated in 
Fig. |4l To observe the sign change experimentally, one 
can add an extra control qubit and only apply erf if the 
control qubit is in the state |0). Explicitly, the recipe is 
as follows, (i) Start from the state 

(|0) + |1)) (E) (|1010) + |0101))/2. (14) 

where the first factor is the state of the control c^ubit and 



(v) Measure the state of the control qubit in the basis 
(±|0) -I- |l))/\/2. Due to the change of sign of the first 
term in step (iii), the control qubit should be found in 
the state (-|0) + \l))/V^ rather than (|0) + |1))/a/2. 

The controlled phase gate required in steps (ii) and 
(iv) can be implemented (up to an overall unimportant 
minus sign) by the following pulse sequence if two Ryd- 
berg levels are available: (i) A tt pulse between |r) and 
the level |c) encoding the state of the control qubit, (ii) a 
27r pulse between |1) and |r'), and (iii) a n pulse between 
|c) and |r). If only one Rydberg level is available, we pro- 
pose to use the pulse sequence illustrated in Fig. \5\ The 
composite pulse (2) in Fig. ^a) consists of three pulses 
with the following parameters, where (f)i is the phase of 
pulse i, fli is the norm of the Rabi frequency of pulse i, 
Ai = uJi — {Er — Ei)/h, uji is the frequency of pulse i, 
{Er — El) is the energy of the Rydberg level relative to 
the energy of level |1), and t is the duration of each of 



the pulses: 

(/.I = 0, 

02 = 7r/2, ^2t 



Aii = 0, 



(18) 



27r sin(\/2r2t) 
1 + cos2(V2m) 
2 V2vr cos 



f^si = nt, 



l + cos2(\/2m) 



A,,i = 0. 



(19) 



(20) 



Let us first consider the case of one atom in |c) and no 
atom in |1). The starting point for the composite pulse 
is then the state with one atom in \r) and no atom in |1), 
which corresponds to the north pole of the Bloch sphere. 
As shown in Fig. [51Jc) , pulse 1 rotates the state an angle 
nt around the x-axis, pulse 2 rotates the state an angle 
^/i^hFfTJA^ = 2tt around the dashed axis in the 
figure, and pulse 3 rotates the state an angle fit around 
the negative a;-axis. As a result, we return to the initial 
state, except that the phase has changed by 



A^io = TT - 



•\/27rcos 



1 + cos2(\/2m) 



(21) 



For the case with no atom in |c) and one atom in |1), 
we start from the south pole but otherwise the picture is 
the same, and the phase of the state is again changed by 
A0O1 = A^io. 

When |c) and |1) are both populated by one atom, we 
start from the north pole, but the Rabi frequency is en- 
hanced by a factor of compared to the values stated 
in Eqs. ([T^ . ([TO]), and (piU) . As a consequence, the sec- 
ond pulse is no longer a 27r pulse. The parameters are, 
however, chosen such that the rotation axis of the second 
pulse (the dashed axis in the figure) passes through the 
point representing the state after the first pulse, and we 
thus again end up on the north pole after the pulse se- 
quence is completed. In this case the phase shift of the 
state evaluates to 



A^n = 



\/27r(l - cos{V2nt)) 
1 + cos2(\/2m) 



(22) 



If |c) and |1) are both unoccupied, the pulses have no 
effect, and the phase change is A6'oo = 0. 

The phase changes obtained through the 
pulses shown in Fig. [Sfb) add to A^oo, 



V27rcos(V2m)/(l + cos2(V2m))i/ 
A9oi, and 2\/27r cos(\/2r2i)/(l -I- cos' 
A6'ii. Choosing fit according to 

cos 



to A6'io and 



to 



■(\/2m))i/2 



(b) |r) 

|e) — \e') 



I -c( \3) |4> 

" ' |0) 




(23) 



FIG. 5. (Color online) Controlled phase gate with a single 
Rydberg level, (a) A tt pulse (1) transfers the population in 
|c) to |r), a composite pulse (2) is applied between |r) and 
which imprints phases but does not alter the populations 
of the levels, and finally a second n pulse (3) transfers the 
population in |r) back to |c). (For simplicity, only the term 
1 1)1 1010) is shown.) (b) Appropriate single qubit phase shifts 
are applied by coupling |c) and |1) off-resonantly to excited 
levels, (c-d) The composite pulse sequence illustrated on the 
Bloch sphere for the case of a total population in jc) and jl) 
of one atom (c) and a total population in |c) and |1) of two 
atoms (d). Note that the Rabi frequency is enhanced by a 
factor of \/2 in the latter case compared to the former case. 
See text for further details. 



the total phase shifts then take the values 0, tt, tt, and tt, 
respectively. 

For the encoding illustrated in Fig. [21 we note that 
it is easier to apply erf operations than af operations, 
since the former do not involve application of the Ryd- 
berg blockade, and it is hence easier to move a flux- 
like quasiparticle around a chargelike quasiparticle. The 
price to pay, however, is the difficulty in preparing the 
ground state of a plaquette, which supports generation 
of a chargelike quasiparticle located at the center of the 
plaquette. 



C. Spinons 

As a final example, we consider spinon excita- 
tions of RVB states. The states under consideration 
are (|0100) - |0010))/V2, (-|0100) + |0001))/V2, and 
(10010) - |0001))/V2, which represent spinons at posi- 
tions 1 and 4, 1 and 3, and 1 and 2, respectively. All 
these states are easy to generate in our system because 
the Rydberg blockade is only needed to move precisely 
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one atom from the reservoir state to one of the states rep- 
resenting qubits. In the basis (|0100), |0010), |0001)), we 
can write the states as (1, -1, 0)^/v^, (-1,0,1)^/72, 
and (0, 1, — 1)'^/-\/2, respectively. 

An interesting property of the states is that they only 
span a two-dimensional space. This space is orthogonal 
to (|0100)-f |0010) + |0001))/v^, which can be experimen- 
tally demonstrated by application of two Raman pulses 
and a measurement. We can, for instance, choose the 
Raman pulses to apply the transformation 



1 

3 







1 

1 











- 



(24) 



in the basis (|0100), |0010), |0001)). This transforms 
our three states (1, -1, 0)^/%/2, (-1, 0, 1)'^/V2, and 
(0,1,-1)^/^2 into (0,-\/3, 1)^/2, (0, VS, 1)^/2, and 
(0,0,-1)-^, respectively. A measurement of the state of 
the qubit represented by level |2) should thus always re- 
sult in the outcome '0', while (1,1, l)"^/\/3 is transformed 
into (1,0,0)^, which populates level |2) with certainty. 
Spinon excitations of the above kind with spin singlet 
pairs of two out of three possible spins (and a random 
state of the uncoupled spin) have recently occurred as 
candidate states for reference frame-free quantum key 
distribution i26|,i2||. 



V. CONCLUSION 

In conclusion, we have proposed to prepare and demon- 
strate basic properties of minimum instances of topolog- 
ical matter in collective excitations of a small atomic en- 
semble. State manipulation is achieved by application 
of the Rydberg blockade mechanism and Raman tran- 
sitions between different atomic levels. In particular, a 
number of transformations that may be difficult to carry 
out on individually selected atoms and pairs of atoms in 
an optical lattice can be achieved in our system simply 
by application of uniform Raman pulses on the entire 
sample. We have provided a general criterion to deter- 
mine whether two states with two excitations each can 
be transformed into each other only by application of Ra- 
man processes, and we have used this result to derive sim- 



ple schemes to prepare and manipulate a number of topo- 
logical states. Experiments on Rydberg blockade have re- 
cently demonstrated controlled two-qubit gates between 
two atoms [1^ [s^] , and we hope that the present work 
will provide inspiration for future experiments demon- 
strating few qubit minimum instances of topological mat- 
ter. Using more internal levels of the atoms and possibly 
more atomic ensembles, our scheme can also be extended 
to generation of topological states of a larger number of 
qubits. As discussed further in [g^], Holmium atoms are 
particularly suitable for this purpose, since they have 128 
hyperfine ground state levels. 

Following the proposal in [glj, the states |5ife2 ■ • • &7v) 
prepared in collective states of the order of a thousand 
atoms can be effectively converted into N well collimated 
wave packets with bi photons in the ith packet. Our 
medium may hence serve as a deterministic source of 
multi-mode field states with a potential of carrying the 
topological quantum correlations over long distances. We 
believe that with the potentially high rate of successful 
generation, the ability to readily create states of higher 
complexity than a single plaquette, and also the possibil- 
ity to either absorb the photon states again in another 
ensemble or interfere them with the deterministic out- 
put from another atomic ensemble, the atomic ensemble 
collective coding is a serious candidate for a quantum re- 
peater for quantum communication protocols with toric 
code qubit protection [Toj . 

Finally, we recall that within the collective encoding 
scheme, we use the Rydberg blockade to restrict the oc- 
cupation of the internal register states to the qubit values 
zero and unity. Higher population of these states can also 
be controlled corresponding to higher spin systems and 
oscillators, for which the plaquette systems offer a wide 
range of interesting topological effects - even in the case 
of coherent states of oscillators, corresponding to clas- 
sical fields in optics j3l| , and to entirely non- interacting 
particles driven by independent Raman transitions in our 
atomic ensembles. 
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